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ABSTRACT 


Crack-growth  resistance  in  transformation-toughened  ceramics  is  studied  by  modeling  the 
region  surrounding  an  advancing  crack  tip  as  a  zone  which  has  undergone  a  uniform  dilatational 
phase  transformation.  This  zone  is  allowed  to  evolve  around  the  advancing  crack  tip  under 
conditions  of  increasing  far-field  load  while  the  tip  is  maintained  at  a  critical  stress  intensity 
necessary  for  fracture.  This  procedure  leads  to  the  surprising  conclusion  that  maximum 
toughening  occurs  for  Finite  amounts  of  crack  advance.  \ 


INTRODUCTION 

^~>The  discovery  of  enhanced  fracture  toughness  in  zirconia-enriched  ceramics  has  led  to  a 
flurry  of  experimental  and  theoretical  analyses.  It  has  been  well  established  that  the  high  stresses 
near  a  crack  tip  can  cause  small  zirconia  particles,  typically  one  micron  or  less  in  diameter,  to 
undergo  a  phase  transformation  from  a  tetragonal  to  a  monoclinic  crystal  structure.  ^  The 
unconstrained  transformation  can  be  decomposed  into  a  four  percent  dilatation  and  a  sixteen 
percent  shear  strain.  However ,  due  to  the  elastic  constraint  of  the  matrix  phase,  the  particles  in  the 
composite  generally  transform  with  twin  bands  of  alternate  character  resulting  in  an  overall 
transformation  strain  considerably  less  than  sixteen  percent.  Consequently,  the  strain 
transformation  of  particles  embedded  in  the  composite  is  usually  assumed  to  be  dilatant. 

The  transformed  region  has  the  remarkable  property  of  allowing  stable  crack  growth  to 
occur  in  the  composite  material,  whereas  the  unreinforced  ceramic  exhibits  no  such  behavior. 
Stable  crack  growth  on  the  order  of  several  millimeters  has  been  observed,  and  R-curves  (i.e. 
stress  intensity  versus  crack  extension)  for  a  variety  of  zirconia-reinforced  composites  have  been 
measured. 
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Theoretical  studies  of  transformation  toughening  including  those  of  McMeeking  and 
Evans  (1982),  Budiansky  et  al.  (1983),  Amazigo  and  Budiansky  (1987),  Rose  (1986a)  and 
Lambropoulos  (1986)  have  considered  only  steady-state  toughening,  wherein  a  semi-infinite  planar 
crack  is  surrounded  by  a  semi-infinite  zone  of  dilatation,  as  shown  in  Fig  la.  Using  a  primitive 
analysis,  Rose  (1986b)  also  examined  the  effect  of  allowing  an  initial  transformation  zone  to 
grow  with  an  advancing  crack  tip.  The  objective  of  this  study  is  to  present  a  complete  analysis  of 
the  growing  crack  by  solving  for  the  initial  transformation  zone  and  studying  its  growth  around  a 
crack  tip  advancing  under  increasing,  far  field,  Mode-I  loading.  This  analysis  will  reveal  the 
surprising  result  that  maximum  toughening  occurs  for  a  finite  amount  of  crack  advance. 

INITIAL  ZONES 

The  modeling  approach  adopted  in  this  study,  consistent  with  previous  work,  is  to  assume 
that  a  transformation  zone  surrounding  the  initial  crack  tip  has  undergone  an  irreversible 
transformation  dilatation  of  strength  c0p,  as  shown  in  Fig.  lb,  where  c  is  the  zirconia  particle 

volume  fraction  and  qJ  is  the  unconstrained  particle  dilatation.  The  crack  will  be  taken  to  be  semi¬ 
infinite  and  planar.  The  transformed  and  untransformed  regions  will  be  assumed  to  have  the  same 
elastic  moduli,  and  to  be  under  plane  strain. 

Since  typical  zone  sizes  are  of  the  order  of  20  pm  or  less,  small  compared  to  specimen 
dimensions,  the  small  scale  zone-size  approximation  will  be  invoked.  The  stresses  at  distances  far 
from  the  zone,  but  small  compared  to  the  size  of  the  body,  will  be  assumed  to  be  dominated  by  the 
classic  Mode  I  K-field .  For  r  -»  ,  the  stresses  are  given  by 

(jap=^=Fop(<|»  (1) 

where  the  Fap(<j>)  are  well  known  trigonometric  functions  and  K  is  the  applied  far-field  stress 
intensity.  Similarly,  in  the  vicinity  of  the  crack  tip,  r  -♦  0,  the  stresses  are  given  by 

<2) 

where  K^p  is  the  stress  intensity  at  the  crack  tip. 
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The  full  transformation  will  be  assumed  to  occur  when  the  mean  stress,  cm  =  0^/3,  attains 
a  critical  value  .  This  corresponds  to  the  "supercritical"  transformation  case  considered  by 
Budiansky  et  al.  (1983).  The  value  of  as  discussed  by  Evans  and  Cannon  (1986),  will  depend 

on  the  stress  and  temperature  history  of  the  composite.  Other  transformation  criteria  have  also 
been  proposed  (Lambropoulos  1986). 

As  shown  by  Budiansky  et  al.,  the  application  of  the  J-integral  to  the  stationa’-y  crack 
results  in  the  conclusion  that  Kyp  =  K.  Accordingly,  if  crack  growth  occurs  at  a  critical  stress 
intensity,  Kup  =  Km ,  where  Km  is  independent  of  the  particle  concentration,  the  initial 
transformation  region  induces  no  toughening.  This  neglects  the  possible  effects  of  other 
toughening  mechanisms  such  as  microcracking. 

To  determine  the  boundary  of  the  initial  transformation  zone  (Fig.  lb),  it  is  necessary  to 
insist  that  the  mean  stress  am  attain  the  value  as  the  boundary  is  approached  from  the  exterior  of 

the  zone.  By  the  superposition  of  stresses,  the  mean  stress  exterior  to  the  zone  boundary  is  the 
sum  of  contributions  from  the  applied  far-field  and  from  dilatations  in  the  transformed  region. 

The  mean  stress  due  to  the  far-field  stresses  of  eqn  (1)  is 


K(l+v) 


(y)  COS(4>/2) 


(3) 


The  mean  stress  contribution  of  the  transformed  zone  can  be  found  by  first  considering  that 
due  to  the  two  small  circular  spots  of  dilatation  shown  in  Fig.  2.  For  spots  of  area  dAo,  strength 
c0p,  and  located  at  Zq  =  XQ+iy0  and  zq  =  x^-iyo,  the  mean  stress  as  found  by  Hutchinson  (1974)  is 


Ec0‘  r  1+v* 


T  1+V1  r.  I 
°m~  ,n-  J  —  v  ^ej 


1 


1 


1811  Ll-vJ-1  ^[VT+V?I 

The  spots  at  z0  and  zq  induce  a  change  in  tip  stress  intensity 

AK“'=drf^>a+s;3'2JdA° 

Rewriting  eqn  (5)  in  terms  of  polar  coordinates  r  and  $,  yields  the  equivalent  form, 


(4) 
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AK,ip  *  3V5T(i-v) cosW2*A°  (6) 

Transformations  in  the  region  <>  <  ti/3  increase  Kiip,  while  those  in  the  region  0  >  tt/3  decrease 
Kup. 

The  effect  of  the  entire  zone  on  the  mean  stress  can  be  calculated  by  integrating  eqn  (4)  over 

the  upper  half  A  of  the  transformed  zone  (Fig.  lb).  The  equation  for  the  zone  boundary  is  then 
obtained  by  adding  the  far  field  mean  stress  to  the  zone  contribution  and  equating  the  sum  to  <4. 

Thus,  the  equation  governing  points  z  =  R^e’*  on  the  boundary  is 

<4  =  cos(<J>/2)  +  J  F(z,z0)dx0dy0  (7) 

where  F(z,zo)  is  given  by  eqn  (4). 

Carrying  out  the  integration  with  respect  to  dxo  and  setting  K  =  Km  leads  to  the  following 
condition  governing  R(<»  when  crack  growth  is  imminent: 

f  L  co  fn  dy0(<j>) 

'HrwJ  ® 

where 

M(<{>,<J> )  =  Re |  - -  +  __  rr-2]!  (9) 

[Vz[Vz  +  ^ZoJ  Vz[Vz  +  ^ojj 


with  z  =  R«>)e1(> ,  Zq  =  R((J>')e,<’  . 
Here  the  parameter  co,  defined  by 


co 


Ec9‘ 
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'  1+v' 
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(10) 


is  a  non-dimensional  measure  of  the  strength  of  the  transformation,  and  the  characteristic  length  L 
is 


L  = 


971 


Km(l+V) 


n2 


'm  J 


Note  that  L  is  the  frontal  intercept  at  0  =  0  of  the  c4  boundary  for  co  =  0. 


(11) 
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The  stress  intensity  at  the  tip  is  equal  to  the  sum  of  the  applied  stress  intensity  and  that 
induced  by  the  presence  of  the  transformation  zone.  Integrating  eqn  (5)  over  upper  half  of  the  zone 
and  adding  the  applied  far-field  intensity  results  in  the  equation 

T 

for  Kup.  Integrating  with  respect  to  dxo  results  in  the  equation 


Ktip  K  2co  f7tfRy1/2 

Kr  =  i£-9FULJ  C0S»,J> 


dy0(4>) 


(13) 


But,  for  the  initial  crack,  K^p  «  K,  and  for  crack  growth  to  occur  KtiP  =  Km .  Hence,  eqn  (13) 
with  Klip/Km=  K/Km=  1  serves  as  a  check  on  the  solution  of  (8)  for  the  initial  boundary  curve 

R(0). 

The  procedure  used  to  solve  eqn  (8)  for  R«>)/L  is  described  in  appendix  A  and  the  results 
for  various  values  of  co  ranging  from  0  to  30  are  shown  in  Fig.  3.  The  effect  of  increased  to  on 
the  relative  shapes  of  the  zones  can  be  seen  in  Fig.  4,  where  for  each  to  all  distances  have  been 
normalized  by  the  frontal  intercept  R(0). 

Figure  5  shows  the  initial  zone-heights  of  Fig.  3,  together  with  the  zone-heights 
corresponding  to  steady-state  crack  growth  calculated  by  Budiansky  and  Amazigo  (1987).  For 
to  =  30,  the  steady  state  zone  height  becomes  infinite,  corresponding  to  the  "lock-up”  -  i.e. 
infinite  toughening  -  discovered  by  Rose  (1986b) 


GROWING  CRACKS 


Once  the  initial  zone  has  been  found,  it  is  possible  to  contemplate  its  growth  around  an 
advancing  crack  tip.  As  a  growing  crack  moves  into  the  body,  material  in  the  vicinity  of  the  tip 
attains  the  critical  mean  stress  and  transforms,  while  due  to  the  irreversibility  of  the  transformation 
a  wake  region  is  left  behind.  Along  a  frontal  portion  of  the  transformed-zone  boundary  the  mean 
stress  criterion  is  satisfied,  while  on  the  wake  portion  of  the  boundary  the  mean  stress  will  have 
dropped  below  om.  According  to  eqn  (6)  material  in  the  transformed  region  which  lies  to  the  left 
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of  the  radial  line  running  through  the  tip  at  the  angle  rc/3  reduces  the  stress  intensity  at  the  tip.  To 
continue  driving  the  crack  forward,  the  applied  stress  intensity  must  be  incremented. 
Consequently  to  solve  the  growing  crack  problem,  both  the  stress  intensity  and  zone  shape  must  be 
found  as  functions  of  crack  extension. 

The  upper  half  of  the  instantaneous  zone  around  a  growing  crack  is  shown  in  Fig.  6a.  The 
boundary'  is  modeled  by  three  segments;  active,  passive,  and  residual.  The  active  segment  AB  is 
that  portion  of  the  boundary  where  the  mean  stress  has  just  reached  o„-  The  mean  stress  on  the 
remainder  of  the  boundary  has  dropped  below  c^.  The  residual  segment  CD  is  the  portion  of  the 

initial  boundary  left  behind  with  the  first  increment  of  crack  growth.  The  intermediate  passive 
portion  BC  is  a  growth  dependent  piece  connecting  the  active  and  residual  segments.  The  passive 
portion  is  comprised  entirely  of  the  end  points  B  of  previous  active  segments,  several  of  which  are 
shown  in  the  sketch. 

The  solution  for  the  growing  crack  problem  involves  adapting  the  three  segment  approach 
to  a  series  of  finite  crack  increments.  In  the  limit  of  continuous  crack  advance,  the  passive  segment 
provides  a  smooth  connection  from  the  residual  portion  to  the  instantaneous  active  segment. 
However  for  a  series  of  finite  crack  increments,  a  piecewise  linear  approximation  to  the  passive 
zone  can  be  constructed  by  connecting  the  residual  segment  to  the  currently  active  segment  with  a 
series  of  straight  lines  running  through  the  end  points  of  previous  active  pieces.  In  the  limit  of 
infinitesimal  crack  increments  the  approximate  boundary  should  coincide  with  the  actual  passive 
segment. 

The  setup  for  a  typical  growth  increment  Aa  is  shown  in  Fig.  6b.  The  active  segment  AB 
is  described  with  respect  to  moving  crack  tip  coordinates  r  and  <}>,  where  $  is  assumed  to  span  the 
angular  interval  from  0  to  an  unknown  angle  a.  For  the  growth  increment  Aa  the  passive  segment 
extension  is  modeled  by  using  a  straight  line  to  connect  the  end  of  the  instantaneous  active 
segment,  B,  to  B’,  the  end  of  the  previous  active  segment.  Tangency  between  the  active  and 
passive  segments  at  B  is  enforced.  For  the  first  crack  increment,  the  passive  segment  connection 
with  tht  initial  zone  shape  must  also  be  found.  As  shown  in  Fig.  6b,  the  passive-residual 
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boundary  at  C  will  be  described  by  the  unknown  angle  P,  measured  with  respect  to  the  initial  crack 
position.  Tangency  between  the  residual  and  passive  segments  will  also  be  enforced  at  point  C. 

The  analysis  presented  for  the  initial  crack  problem  can  be  applied  to  the  growing  crack 
configuration  with  some  slight  modifications.  The  initial  zone  will  be  allowed  to  grow  with  an 
advancing  crack  tip  under  the  dual  requirements  of  satisfying  the  critical  mean  stress  condition  on 
the  active  zone  boundary  and  maintaining  the  tip  stress  intensity  at  Km. 

The  mean  stress  criterion,  eqn  (7),  is  enforced  with  the  understanding  that  the  angle  6  is 
restricted  to  the  interval  0  <  <f>  <  a.  Introducing  the  nondimensional  toughening  parameter 
A  =  K/Km,  carrying  out  the  integration  with  respect  to  dxo,  and  regrouping  eqn  (7)  results  in  the 
expression 

to  f"  dy0«>') 

1  =  j  cos($/2)  -  ^  Jo  M «>,<>')  (14) 

governing  the  active  zone  boundary  R(<»/L.  All  distances  in  eqn  (14)  are  measured  from  the  crack 
tip  and  the  integration  extends  over  the  entire  boundary  0  <  <  7t . 

The  tip  stress-intensity-factor  can  be  maintained  at  Km  by  enforcing  eqn  (13),  which  now 

requires 

.  ,  2ci)  r’VRV1'2  dyo(« 

A=1+97J0IlJ  cos  m)  —  (,5) 

The  system  (14)  and  (15)  constitute  a  nonlinear  integral  equation  and  a  scalar  equation  for 
R«>)/L  and  A.  For  the  initial  crack  increment,  (14),  (15)  and  tangency  conditions  at  a  and  p,  can 
be  solved  for  R(<}>)/L  in  (0  <  <>  <  a),  A,  a  and  p.  For  a  series  of  subsequent  crack  increments 
R(d)/L,  A,  and  a  are  found  repeatedly  to  generate  the  zone  shape  and  the  crack-growth  resistance. 
The  solution  procedure  for  each  crack  increment  is  outlined  in  Appendix  B. 

Growing  zone  shapes  for  co  =  5  and  10  are  shown  in  Fig.  7.  The  dashed  curves  indicate 
the  positions  of  active  segments  for  various  amounts  of  crack  extension;  the  innermost  curve  is  the 
initial  boundary.  The  R-curves,  plots  of  A  versus  Aa/L,  for  o>  =  5,10  along  with  their  respective 
steady  state  toughening  asymptotes,  as  found  by  Budiansky  and  Amazigo  (1987),  are  shown  in 
Fig.  8.  The  results  of  Figs.  7  and  8  are  quite  unexpected.  The  zone  height  ,  H/L,  and  the 
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toughening,  A,  both  overshoot  their  steady-state  levels  for  finite  amounts  of  crack  advance  before 
approaching  them  asymptotically  from  above. 

The  R-curves  of  Fig.  8  are  qualitatively  consistent  with  some  available  experimental 
measurements.  A  number  of  investigators,  including  Swain  (1983),  and  Swain  and  Hannink 
(1984),  have  reported  R-curves  which  exhibit  peaks  in  toughness.  However,  the  corresponding 
peak  in  zone  height  has  not  been  reported.  Even  though  the  actual  zone  boundary  between 
transformed  and  untransformed  regions  occurs  over  a  diffuse  region,  the  issue  of  zone  widening 
should  be  explored  experimentally. 

R-curves  and  transformation-zones  have  been  calculated  for  various  values  of  co,  and  the 
results  for  the  peak  toughening  Ap  and  peak  zone  height  Hp  are  shown  in  Table  1.  Also  shown  are 
the  crack  extensions  Aa(Ap)  and  Aa(Hp)  necessary  to  reach  each  peak.  Both  Ap  and  Hp  are  seen  to 
increase  dramatically  in  the  interval  19  <  co  <  20.  Denoting  the  steady  state  toughening  ratio 
Ks/Kn„  as  calculated  by  Budiansky  and  Amazigo  (1987),  by  As,  we  compare  the  toughening 
predictions  for  the  steady  state  and  growing  crack  calculations  in  Fig.  9,  which  show's  (A*)-1  and 
(Ap)’1  versus  co.  The  implication  for  a  growing  crack  is  "lockup"  for  co  greater  than  20.4  .  The 
steady  state  configuration  (Fig.  la)  can  not  be  reached  from  the  growth  of  an  initial  zone  for  co  > 
20.4  . 

It  may  be  noted  from  Fig.  8  that  a  tiny  secondary  maximum  appears  in  the  R-curve  for 
co  =  10  .  For  values  of  co  greater  than  10  (but  less  than  20.4)  additional  peaks  of  decreasing 
amplitude  were  found  to  appear  in  the  R-curves  as  the  toughness  decayed  in  an  oscillatory  fashion 
to  its  steady-state  magnitude. 

By  considering  the  initial  slope  of  the  R-curves,  it  is  possible  to  define  an  initial  tearing 
resistance  parameter 


d(K/Km)  E  ( l+v'i  dJ 


Here  t  represents  the  rate  of  toughening  increase  as  the  crack  begins  its  initial  growth  and  is 
plotted  versus  co  in  Fig.  10.  A  comparison  can  be  made  with  the  approximate  results  of 
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Hutchinson  (1986)  who  neglected  the  effect  of  the  transformation  on  the  zone  boundary  in  his 
study  of  initial  tearing  resistance.  Figure  1 1  shows  a  comparison  of  both  sets  of  calculations  for 
small  co.  In  the  limit  co  -» 0,  the  predictions  coincide. 

CONCLUDING  REMARKS 

In  the  presence  of  transforming  panicles,  the  resistance  to  crack  growth,  as  measured  by 
the  applied  stress  intensity  K,  reaches  a  maximum  at  a  Finite  amount  of  crack  growth.  This 
maximum  exceeds  the  steady-state  toughness  by  an  amount  that  varies  with  the  transformation 
intensity  parameter  co.  Furthermore,  "lock-up"  (i.e.  infinite  toughening),  occurs  at  a  critical  value 
of  co  =  20.4,  which  is  substantially  smaller  than  the  value  co  =  30  corresponding  to  steady-state 
toughening.  The  potential  benefits  of  transformation  toughening  must  be  assessed  on  the  basis  of 
transient  crack  growth  since  the  steady-state  condition  seriously  underestimates  the  toughening 
effect  of  transforming  panicles. 
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APPENDIX  A 

The  initial  zone  boundary,  R«>),  can  be  found  by  solving  the  nonlinear  integral  equation 
(8).  However  before  proceeding,  it  is  instructive  to  examine  the  solution  for  co  =  0  given  by  eqn 
(3)  and  shown  in  Fig.  3: 

R(0)  =  L[l/2  +  1/2  cos($)]  (Al) 

Due  to  Mode-I  symmetry,  the  boundary  at  4>  =  0  intersects  the  axis  ahead  of  the  tip  normally. 
However,  at  <$>  =  n,  the  boundary  runs  through  the  crack  tip  tangentially  to  the  crack  face. 

For  co  *  0,  the  possibility  that  the  boundary  at  4>  =  7t  detaches  from  the  crack  tip  must  be 
admitted.  An  analysis  of  the  integral  in  eqn  (8)  shows  that  if  R(tt)  *  0,  the  boundary  intersects  the 
crack  faces  normally.  An  expansion  for  R(4>)  meeting  the  boundary  requirements  at  0  and  rt  is 

N 

R«0  =  X  ®n  cos(n<»  (A2) 

n=0 

Substitution  of  equation  (A2)  allows  eqn  (8)  to  be  rewritten  in  the  form, 

G(ao,a1,...,aN,0)  -1=0  tA3) 

where  the  function  G  is  the  right  hand  side  of  eqn  (8). 

In  this  study,  the  solution  for  the  (N+l)  unknown  coefficients  was  accomplished  by 
collocating  eqn  (A3)  at  N+l  equally  spaced  points,  4>  =  7tn/(N+l)  where  n  =  (0,1,2,...,N),  in  the 
interval  0  <  <>  <  7t .  A  Newton-Raphson  iterative  technique  was  then  used  to  solve  for  the  N+l 
unknown  coefficients.  The  integral  in  eqn  (A3)  was  evaluated  by  Gaussian  quadrature. 
Convergence  was  assumed  when  the  relative  change  in  successive  iterations  of  each  of  the 
unknowns  was  less  than  0.001.  For  the  initial  zones  shown  in  Fig.  3,  it  was  found  that  a  10  term 
expansion  was  sufficient  to  obtain  a  highly  accurate  solution.  As  an  additional  check  on  the 
accuracy  of  the  overall  shape,  the  constraint  on  R($)/L,  eqn  (13)  with  Kup/K  =  1,  was  evaluated 
and  found  to  hold  to  within  H)-5  for  each  of  the  zones  of  Fig.  3. 
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APPENDIX  B 


The  solution  for  the  unknowns  R(<{>)»  a,  and  A  after  a  typical  crack  increment  (Fig-  6b) 
involves  solving  the  integral  equation  (14),  the  scalar  equation  (15),  and  meeting  the  tangency 
condition 

^  (®  )  I  passive  =  ^  (®)  I  active  (B  1 ) 


For  the  initial  crack  increment  the  system  must  be  supplemented  with  the  unknown  P  and  the 
additional  tangency  equation 


fwl 


_  dy  | 

passive  ~  ^  v'-'/  ■  residual 


(B2) 


An  appropriate  expansion  for  the  radius  of  the  instantaneous  active  segment  which  meets 
Mode-I  symmetry  conditions  is 


N 

(B3) 

0*0 

where  the  T2D  are  the  even  Tchebyshev  polynomials.  Substituting  this  expression  into  (14),  (15), 
and  (Bl)  results  in  a  set  of  equations  for  the  unknowns  (ao,ai,...,a>sj),  a, and  A. 

The  solution  procedure  will  be  described  for  the  general  growth  increment  as  shown  in  Fig. 
6b.  For  the  initial  increment,  the  process  is  the  same  with  P  and  (B2)  added  to  the  list  of 
unknowns  and  equations.  For  a  given  Aa,  a  system  of  N+3  equations  was  generated  by  enforcing 
equations  (15),  (Bl),  and  collocating  eqn  (14)  at  N+l  equally  spaced  points  <|>  =  na/(N+l) 
(n=0,l,2,...,N+l).  The  equations  were  written  in  residual  form  and  Newton-Raphson  method 
was  used  to  find  the  solution.  Gaussian  quadrature  was  used  to  evaluate  the  integrals  in  eqns  (14) 
and  (15).  Convergence  to  a  solution  was  specified  by  a  relative  change  of  less  than  0.001  between 
iterations.  A  6  term  expansion  for  R(4>)  was  found  sufficient  to  obtain  a  solution  at  all  stages  of 
growth. 

Initial  growth  increments  had  to  be  small,  <  0.05/L,  to  capture  the  early  zone  shape.  As 
crack  extension  proceeded,  it  was  found  that  the  increment  size  could  be  increased.  The  number  of 
increments  necessary  to  grow  the  crack  to  peaks  in  height  and  toughening  varied  greatly.  For 
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small  values  of  to  less  than  50  increments  were  sufficient.  However  for  values  of  to  between  17 
and  20,  several  hundred  increments  were  necessary  to  reach  Hp  and  Ap.  During  the  growth 
process,  the  mean  stress  exterior  to  the  zone  was  checked  to  confirm  that  it  was  less  than  <£■  The 

final  results  shown  herein  were  checked  by  calculations  with  the  increment  steps  halved. 
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TABLE  1 


(0 

Ap 

Hp/L 

Aa(Ap) 

Aa(Hp) 

5 

1.29 

1.03 

5.5 

2.4 

10 

1.80 

1.91 

10.25 

5.6 

15 

3.07 

5.17 

29.4 

18.4 

17.5 

5.06 

13.6 

81.0 

55.0 

19.5 

16.3 

351 

2400 

1875 

20.0 

74.3 

2764 

20,200 

15,900 
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Fig.  la  Steady-state  transformation  zone  surrounding  a  crack  tip 


Fig.  lb  Initial  transformation  zone  around  a  stationary  crack 
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Fig.  2  Symmetrically  located  spots  of  dilatation 


Fig.  6b  Zone  extension  model  used  to  find  the  instantaneous  active  zone  and  passive 
segment  growth 
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for  co  =  5,10 


